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Abstract. We consider the energy spectrum of cosmic-rays (CRs) from a purely phenomenological point of view and investigate 
the possibility that they all be produced by the same type of sources with a single power-law spectrum, in E~\ from thermal 
to ultra-high energies. We show that the relative fluxes of the Galactic (GCR) and extra-galactic (EGCR) components are 
compatible with such a holistic model, provided that the index of the source spectrum be x m 2.23 ± 0.07. This is compatible 
with the best-fit indices for both GCRs and EGCRs, assuming that their source composition is the same, which is indeed the 
case in a holistic model. It is also compatible with theoretical expectations for particle acceleration at relativistic shocks. 
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1. Introduction 

Despite considerable observational and theoretical progress 
over the last decades, the origin of cosmic-rays (CRs) re- 
mains uncertain at all energies and a subject of intense de- 
bate. The very existence of ultra-high energy CRs (UHECRs) 
challenges the most efficient particle acceleration mechanisms 
thought to be at work in powerful astrophysical sources, and 
the CR source model often assumed for low energy particles 
(involving diffusive shock acceleration in isolated supernova 
remnants) suffers from a number of persistent problems that 
shed some doubts about its general validity (e.g. Parizot, et al., 
2001, Parizot 2005). The only well-established facts about low- 
energy CRs are derived from the detailed and joint study of 
their composition and spectrum between, say, 100 MeV/n and 
100 GeV/n, in the light of CR propagation models involving 
energy losses, in-flight spallation reaction and charged parti- 
cle diffusion and confinement in a magnetised extended halo 
(e.g. Strong and Moskalenko, 2001). They can be summarised 
as follows: phenomenologically, i) the source spectrum is con- 
sistent with a power-law, E~ x , of index x - 2.2-2 A (e.g. Jones, 
et al., 2001; Ptuskin et al., 2005), ii) the source composition 
is compatible with the acceleration of the average interstellar 
medium, somewhat modified by specific selection mechanisms 
and local enrichment from massive star's ejecta (e.g. Meyer et 
al., 1997; Ellison, et al., 1997), iii) the diffusion coefficient of 
the relativistic nuclei increases with their energy-to-charge ra- 
tio as (E/Zf, where /? * 0.3-0.6 (e.g. Jones, et al., 2001), and 
iv) the height of the confining halo is relatively large, of the or- 
der of 5 kpc above and below the Galactic plane (Moskalenko, 
et al., 2001). As for the UHECRs, one can typically distin- 
guish between two types of phenomenological scenarios (for 
astrophysical scenarios with roughly homogeneous extragalac- 
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tic sources): either pure proton sources with a steep source 
spectrum in E~ 26 or E~ 2J or sources with a mixed composi- 
tion (similar to the low-energy cosmic-rays, say) and a harder 
source spectrum, in £~ 2 2 or E~ 2 3 , typically (see Allard et al., 
2005a,b, and refs. therein). 

It is customary to distinguish between several regions of 
the CR spectrum. Some of these distinctions are artificial or re- 
lated to the history of CR science or to the specific detection 
technique used in different energy ranges, while others have 
more phenomenological or theoretical grounds, in relation with 
observed or expected structures in the spectrum (e.g. due to 
solar modulation, spatial diffusion, contributions of individual 
sources, loss of confinement, energy losses...). However, de- 
spite the spectral features known as the knee(s) and ankle, the 
most striking property of the CR spectrum is its overall regu- 
larity and coherence over (at least) 12 orders of magnitude in 
energy and 32 orders of magnitude in differential flux. This is 
quite extraordinary and unique for a non-thermal phenomenon. 
It invites one to consider cosmic rays as a global, unified phe- 
nomenon taking place in the universe, and investigate how far 
it is possible to follow such a road. 

However, even within a unified point of view, a distinction 
cannot be avoided between CRs of Galactic and extragalactic 
origin. Indeed, low-energy cosmic rays (at least up to Te V ener- 
gies) are known to be accelerated in our Galaxy, and the lower 
gamma-ray emissivity from neutral pion decay observed in the 
large Magellanic cloud (LMC), as compared to the Milky Way, 
implies that the CR density inside our galaxy is larger than in 
the LMC (Sreekumar, et al., 1993). This would not be the case 
if the responsible CRs had an extragalactic origin. At the other 
end of the spectrum, on the contrary, CRs with an energy larger 
than, say, 10 19 eV have too large gyroradii to be confined by the 
Galactic magnetic fields, and their observed (rough) isotropy 
thus implies an origin outside (of the disk) of the Galaxy. At 
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least two components of CRs are thus observed, one Galactic 
and one extragalactic, with a transition probably around the so- 
called ankle, at ~ 3 10 18 eV (e.g. Allard, et al., 2005a,b) or 
below (Berezinsky et al., 2004). 

However, the fact that the low-energy and high-energy CRs 
that we observe have a different origin in space does not imply 
that they have a different origin in nature. Following a unifi- 
cation perspective looking for as few astrophysical processes 
as possible behind the CR phenomenon, we can investigate the 
possibility that all CRs be actually produced by a single type 
of sources, from thermal energies up to a maximum energy, 
Emax > 10 20 eV, that still remains to be determined. We refer to 
such a model as a holistic model. Since, as recalled above, the 
CR sources are quite uncertain at all energies, this opens the 
possibility that one single model may solve both problems. 

Different versions of holistic models have been proposed 
before, notably with gamma-ray bursts as the source of all CRs, 
within the cannonball model (Dar and de Rujula, 2004). In a 
recent version (Dar, 2004), CRs below and above the knee are 
accelerated by the same astrophysical objects, but not by the 
same mechanism nor with the same spectrum. In another ap- 
proach (Dar and Plaga, 1999), all the observed CRs originate 
from our galaxy and its halo, with a negligible contribution of 
other galaxies, even at ultra-high energies. In this Letter, we 
consider a different (strong) version of holistic models, where 
not only the same source, but also the same acceleration mech- 
anism is operating over the whole spectrum of CRs, with the 
simplest source spectrum, namely a single power-law, in E~ x . 
This is a purely phenomenological study, where the actual CR 
sources and acceleration mechanism are left unspecified. We 
simply assume that some mechanism produces a CR power-law 
spectrum in some sources located inside galaxies, and investi- 
gate the viability of such a generic model in view of the avail- 
able CR data. We show that the resulting relative normalisation 
of Galactic and extragalactic components are as observed for a 
very natural value of the unique free parameter of the model, 
namely the logarithmic slope of the power-law, x. 

2. Relative normalization of Galactic and 
extragalactic CR components 

In order to compare the observed relative normalisation of 
the galactic and extragalactic CRs components with the one 
predicted by the holistic model, we choose a reference en- 
ergy in each component and compute the corresponding CR 
density, respectively inside our Galaxy and in the whole uni- 
verse. We choose Eq = 1 GeV as the reference value for 
(low-energy) Galactic cosmic -rays (GCRs), because CR phe- 
nomenology (including propagation effects) is well known and 
constrained by composition measurements around that energy. 
The CRs injected in the Galaxy are resonantly scattered by 
magnetic field inhomogeneities and remain confined for a time 
Tconf (E) that depends on their energy. At very high energy, say 
at E £ 5 10 18 eV, the CR gyroradius is larger than the magnetic 
field coherence length and even becomes comparable to the 
size of the confining region, so that the CRs injected inside the 
Galaxy leak out quickly and diffuse away in the whole universe. 
Thus, as recalled above, most of the high-energy CRs observed 



on Earth have an extragalactic origin. We choose E\ = 10 19 eV 
as the reference value for extragalactic cosmic-rays (EGCRs), 
because: i) it is high enough for the corresponding CRs to have 
a diffusion radius (over the energy loss time scale; i.e. a mag- 
netic horizon, see Deligny et al. 2004, Parizot 2004) larger 
than the distance between typical galaxies, so that their den- 
sity should be roughly uniform in the universe, and ii) it is low 
enough for a reliable estimate of the corresponding CR density 
to be obtained in a simple way, both observationally (since the 
flux is not yet dominated by statistical uncertainties) and theo- 
retically (since E\ is below the GZK regime where stochastic 
propagation effects and the local source distribution have a no- 
ticeable influence on the observed fluxes). 

2.1. GCR density 

Let us assume that CRs are produced in our Galaxy with an 
average injection rate, in s -1 GeV -1 , 



Q(E) = Q 



E 



(1) 



where the normalisation, go, relates to the global source power. 
In practice, it does not matter where and when the CRs are 
injected in the Galaxy, as long as the granularity of the sources 
in space and time is small compared to the diffusion radius and 
confinement time of the CRs in the Galaxy. This amounts to 
assuming that the observed CRs are not significantly different 
from what would be observed at another time and somewhere 
else in the Galaxy, as usually assumed in CR studies. 

At low-energy (and in a steady state), the local density of 
the GCRs is then simply given by the total number of CRs in- 
jected during the confinement time, T conf , divided by the con- 
finement volume, V con f : 



, , E \ X r conf (E) 



(2) 



2.2. EGCR density 

At high-energy, many (extragalactic) sources contribute to the 
local density of EGCRs, because of the dilution of the particles 
over large volumes in the universe. The dilution volume de- 
pends on the transport properties in the extragalactic medium. 
In the absence of magnetic fields, it is merely a sphere of radius 
ct, where r is the "lifetime" of the EGCRs (if shorter than the 
age of the source), while in a diffusive transport regime with 
diffusion coefficient D, the diffusion radius is given by V6Z)t. 
For an order of magnitude calculation, one may consider that 
the particles injected at an energy E "disappear" when they 
have lost a substantial fraction of their energy, and thus choose 
the EGCR lifetime as the energy loss time, Ti 0SS (E), defined by 
the energy drift equation dE/dt = E/t\ 0SS (E). At the reference 
energy, ti oss (£i) 4.4 10 9 years, which is large enough to en- 
sure an overlap of the dilution spheres around the galaxies. 

Then, considering that the EGCRs accumulate uniformly 
throughout the universe for a time ti oss , the density of EGCRs 
is obtained as the number of CRs injected by a galaxy dur- 
ing that time, multiplied by the density of galaxies in today's 
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universe, n ga i. In practice, galaxies are of different types and 
one may expect that they inject CRs (roughly) proportionally 
to their star formation rate. Their contribution to the over- 
all CR density should then be scaled accordingly. Here, we 
take the Milky Way as a reference and thus use the density 
of "Milky Way-like galaxies" in the universe, estimated from 
the ratio of the star formation rate density in the local uni- 
verse, 1.4 lO^Moyr'Mpc- 3 (Wyder et al., 2005) and the 
Galactic star formation rate, 3M yr~' (e.g. McKee, 1989): 
n g ai - 5 10~ 3 Mpc~ 3 . We thus have: 



n EG (E) — <2o -pr 



x riossCE) X n ga i , 



(3) 



to be compared with Eq. 0. 

More rigourously, the EGCR density at energy E is ob- 
tained from the sum of all particles injected at earlier times 
(or redshifts z) at a higher energy such that they have lost just 
enough energy to be observed today at energy E: 



n E c(E) = n ga i X 



dE 

Q(E in (E,t lb ),t lb )—^-dt lh , 
dE 



(4) 
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Fig. 1. Cosmological factor, I(E\,x), defined in Eq. (|8), as a 
function of the CR source spectrum index, x, for £1 = 10 EeV 
and two different choices of the source power evolution, f(z). 



where t\ b is the look-back time, from today (%> = 0) back to 
the onset of EGCR sources {t\ b = T), Ei n (E, t\ b ) is the injec- 
tion energy, at t\ b , of EGCRs observed today at energy E, and 
Q(E, t) = Qo(E/Eo)~ x x /(f), allowing for a possible evolu- 
tion of the source power (/(0) = 1). The injection energy is 
computed by solving the energy loss equation dE/dt = E(E, ?), 
where the latter function includes both expansion losses, 



dE 
~dt 



exp 



E dt 

x — » 

I + z dz 



(5) 



and losses due to the EGCR interaction with the cosmological 
background radiation (e.g. Berezinsky and Grigorieva, 1988): 



dE 
~dt 



E(l+z) 



CMB 



Jo 



((tk)[xE{\ +z)] xln(l -e-P x )dx, (6) 



where (ctk) is the product of the interaction cross section (for 
e + -e~ and pion production) by the corresponding inelasticity, 
the constants K and /? are given by K — {k^Tom 2 c A )l{2nh i c 2 E) 
and /} = m p c 2 /2kftTo, and To = 2.73 K. 
Equation (0} can be rewritten as 



x Tio SS (£) x n gal x I(E, x), (7) 



«eg(£) - 2o I — 



where ti oss = E/E(E,0) is computed at z — and the integral 
I(E, x) gathers all the modifications to the simplified Eq. Q 
due to the cosmological evolution of the sources and the time- 
dependent energy losses: 



I(E, x) = 



E(E,Q) 



r°/E- m \- x dE ia dt 

J„ (f ) -di f(z) dz dz - (8) 



Note that I(E, x) would be equal to \/{x - 1) if E(E, t) were 
constant and the sources did not evolve, i.e. f(z) - 1. If Ho 
is the current Hubble "constant", Qm is the usual normalised 



matter density, Qa is the reduced cosmological constant and 
Qk = 1 - Qm - (the "curvature density"), we have: 



df _ 1 

dl ~ Ho(\+z) 2 



n M (l +Z) + + 



(1 +z) 2 



-1/2 



(9) 



This allows us to compute I(E, x), as shown on Fig. [2 (we 
assume fi M = 0.27, Q A = 0.73 and H - 75 km/s/Mpc). It 
is convenient here to approximate I(E, x) by an exponential 
function of x. To within a few percent in the range of x of inter- 
est, we found I{E\,x) =t 1.4 e -a53 * in the case when f(z) = 1 
(no source evolution) and I{E\, x) ^ 3.7 e~ 014x with a source 
evolution in f(z) = ( 1 + z) 4 (by analogy with the star formation 
rate). 

2.3. GCR/EGCR normalisation 

The predicted value of nc(£b)/«EG(£i) follows straightfor- 
wardly from Eqs. and 0: 



n G (Eo) 



TconfCEo) 
X ■ — - X 



1 



TlossCEl) « ga l Vconf (E )I(Ei , x) 



(10) 



Identifying the left hand side with the measured CR flux ratio 
and using writing I(E\ , x) Ae~ Bx (see above), one can easily 
solve the above equation to determine the value of the only 
free parameter in the model, namely the logarithmic index of 
the (holistic) injection spectrum, x: 



In 



A ngalVconftEo) 



r loss (£,) €> C r(£o) 



Tconf (Eo) <t>CR(£l) 



/in 



Eo 



(11) 



The required value of x is thus given here as a function 
of parameters that can be measured or derived directly from 
astrophysical data. Numerically, we gave above the value of 
TiossCEi) and an estimate of n ga t. The confinement time at 1 GeV 
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if roughly T conf (E ) * 2.4 10 7 years (Connel, 1998). As re- 
called above, recent propagation studies favour halo heights 
of the order of 5 kpc (above and below the Galactic plane). 
Assuming a CR disk radius of 20 kpc, one may thus esti- 
mate the confinement volume at 1 GeV as Vconf(^o) — 71 x 
(20 kpc) 2 x 10 kpc a 1 .3 10~ 5 Mpc 3 . For the measured CR flux 
at Eo we take the differential flux deconvoluted from solar mod- 
ulation (Webber, 1998): <D CR (1 GeV) 0.5 cnrV-^GeV 1 , 
and at E\ we take an average of AGASA (Teshima, et al., 
2004) and HiRes (Zech, et al., 2004) values: <D CR (10 19 eV) a 
210- 28 cirT 2 sr- 1 s- 1 GeV- 1 . Therefore, <D cr (£o)/3>cr(£i) - 
2.5 10 27 . Inserting these values into Eq. (lilt , and assuming ei- 
ther f(z) = (1 + z) 4 , so that A = 3.7 and B - 0.74, or f(z) = 1, 
so that A = 1.4 and B = 0.53, we find: 

x a 2.23 ± 0.07, for/(z) = (1 + z) 4 , 

J (12) 
x a 2.21 ±0.07, for /(z) = 1, 

where the "error bars" account for an indicative uncertainty of 
a factor of 5 on the quantity «g a iV CO nf(£ , o)(Ti oss (£ , i)/T conf (£'o)), 
as well as a possibly different source evolution function. 

3. Discussion 

The result shown in Eq. d!2l > reads as follows: holistic CR 
source models (where CRs of all energies are produced by the 
same sources with a single power-law spectrum) are indeed 
possible, and the corresponding spectral index is determined 
unequivocally from measured parameters: x a 2.2-2.3. This 
is remarkable in several respects. Not only could a sensible 
value of x be found, but the solution is particularly interesting 
from the phenomenological and theoretical points of view. This 
slope is indeed in keeping with common expectations, and in 
remarkable agreement with relativistic shock acceleration pre- 
dictions (e.g. Kirk et al., 2000). It is also right in the range 
allowed for low-energy CRs, as derived from detailed studies 
of the secondary-to-primary composition ratios (cf. Sect.^. 

The above solution is also particularly interesting in rela- 
tion with EGCRs. While the best fit of the highest-energy data 
is usually obtained with a source spectrum in E~ 2 6 (De Marco, 
et al., 2003), this result only holds for pure proton sources. 
Most (if not all) astrophysical sources of CRs would however 
accelerate heavier nuclei just as well, since electromagnetic 
processes only depend on the charged particle rigidity. Quite 
remarkably, it is found that when assuming a similar composi- 
tion for EGCRs and GCRs (as would of course naturally be the 
case in a holistic model), the high-energy data are best repro- 
duced by a power-law source spectrum in E~ 2 2 or E~ 23 , de- 
pending on source evolution (Allard et al., 2005a,b). Thus, the 
source spectrum that makes holistic models viable is precisely 
the one that is consistently favoured at low energy (by GCR 
studies) and, independently, at high energy (by EGCR studies). 

Even though it is simple, the above analysis is quite robust. 
The most uncertain parameters are the CR confinement volume 
at 1 GeV and the density of Milky-Way-equivalent galaxies in 
today's universe. However, even large variations of these pa- 
rameters cannot change the value of x substantially, as a re- 
sult of the huge lever arm between low-energy GCRs and high- 
energy EGCRs. It should also be noted that the determination 



of x does not depend on any assumption concerning the energy 
dependence of the CR confinement time. We simply used the 
value derived from the CR data at Eq = 1 GeV. Interestingly, 
however, it then leads to a determination of the confinement 
time that would allow one to reproduce the observed GCR 
spectrum, in ~ E~ 2J1 : T coni (E) cc E s , with 0.41 < 5 ^ 0.55, 
which is also the range expected from CR diffusion theory. 
More generally, all the phenomenological studies concerning 
GCRs and EGCRs turn out to be equally valid within the holis- 
tic models, because the required spectral index is the one that 
was derived (or needed) at both low and high energies anyway. 

In conclusion, we found that a holistic CR source model 
with a single power-law spectrum in E~ x and x a 2.23 + 0.07 
could account for the CRs at all energies, in keeping with all 
known results about CR phenomenology, and in agreement 
with the main theoretical results for CR acceleration and trans- 
port. However, it should be clear that we did not propose a 
well-defined physical or astrophysical model. We simply made 
a general phenomenological remark that may motivate or en- 
courage further studies of the origin of CRs considered glob- 
ally, as a general phenomenon possibly involving a single pro- 
cess at work over their whole energy range. 

Finally, one can propose another reading of the result ob- 
tained here. Reversing the point of view, one may state that, 
whatever the sources of GCRs, their spectrum actually ex- 
tends up to ~ 10 20 eV, then the contribution of the corre- 
sponding high energy particles escaping from all galaxies is 
just enough to explain the EGCR fluxes and properties (i.e. no 
other sources are needed), provided that the source spectral in- 
dex is ~ 2.2-2.3. Likewise - and perhaps more strikingly-, 
even though we do not know what the EGCR sources are, it is 
very natural to expect that they do not produce only UHECRs, 
and that the latter are just the high-energy end of a continu- 
ous spectrum starting at thermal energies in the sources. Now, 
according to our result, the EGCR source spectrum extends 
down to non-relativistic energies indeed, then the contribution 
of the corresponding low-energy CRs confined within the par- 
ent galaxies is just enough to explain the GCR flux and prop- 
erties as well (i.e. no other sources are needed), provided that 
the source spectral index is ~ 2.2-2.3, which it is inferred to be 
anyway, from (mixed-composition) EGCR studies. 

Acknowledgements. EP thanks Arnon Dar and Paul Sommers for in- 
teresting discussions about the model. 
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